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Abstract 

Using quantum field theory and bosonization, we determine the quantum 
phase diagram of the one-dimensional Hubbard model with bond-charge 
interaction X in addition to the usual Coulomb repulsion U at half-filling, 
for small values of the interactions. We show that it is essential to take 
into account formally irrelevant terms of order X. They generate relevant 
terms proportional to X 2 in the flow of the renormalization group (RG). 
These terms are calculated using operator product expansions. The model 
shows three phases separated by a charge transition at U = U c and a spin 
transition at U = U s > U c . For U < U c singlet superconducting correlations 
dominate, while for U > U s , the system is in the spin-density wave phase 
as in the usual Hubbard model. For intermediate values U c < U < U s , the 
system is in a spontaneously dimerized bond-ordered wave phase, which is 
absent in the ordinary Hubbard model with X = 0. We obtain that the 
charge transition remains at U c — for 1^0. Solving the RG equations 
for the spin sector, we provide an analytical expression for U S (X). The 
results, with only one adjustable parameter, are in excellent agreement with 
numerical ones for X < t/2 where t is the hopping. 
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1. Introduction 



1.1. The model 

The Hubbard model, with nearest-neighbor hopping t and on-site repul- 
sion U has been widely used to study the effects of correlations, as a simpli- 
fied model to describe compounds of transition metals and other systems. 
However, one expects that in any system, in general the hopping between 
two sites depend on the occupation of these two sites, which leads to the 
presence of bond-charge interactions (also called correlated hopping terms) 



in the Hamiltonian [such as X in Eq. (jTJ and /^-v (cj^v +H.c. 



in,- 



T nj. 



For example, in simple systems with one relevant orbital per site, one would 
expect that when electrons are added to one site, the screening of the core 
charge increases, and as a consequence, the wave function of the orbital 
expands and the hopping to the nearest sites increases. In general, any one- 
band effective model derived from more complex Hamiltonians to describe 
the low-energy physics of some system, contains bond-charge interactions. 
In fact, the generalized Hubbard model with correlated hopping terms has 
been derived and used to describe the low-energy physics of intermediate 
valence systems [l|, organic compounds @, IE 0, 11, @, 0]) a Hubbard model 
including lattice vibrations (Hi, c uprate superconductors 0,0, 111! ], and 



more recently optical lattices [12|, [13J, [14[ . This is particularly interesting 



because the parameters can be tuned experimentally in a wide range 1 4J. 1 1 5 



First-principles calculations in transition-metal complexes suggest that the 
correlated hopping terms can be large 16[ . 

As we shall see, the presence of bond-charge interaction leads to quali- 
tatively new physics. One example is that in two dimensions d-wave pair- 

3. 



arc 



18 



ing correlations, which are already present in the Hubbard model 
strongly enhanced in the generalized Hubbard model for the cuprates 
and one obtains ci-wave superconductivity already at the mean-field level 
19]. In one dimension (ID), field-theoretical 2(J 21] and numerical 21, 22 



results show the presence of a spontaneously dimerized bond-ordering wave 
(BOW) and a phase with dominant triplet superconducting correlations at 
large distances that are absent in the ordinary Hubbard model. For special 
values of the parameters, the model with two- and three-body interactions 
has been solved exactly by the Bethe ansatz 23 . 
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The simplest model with bond-charge interaction has been proposed by 



24, 25 . The 



Hirsch motivated by his theory of hole superconductivity 
Hamiltonian can be written as 



H = -t ^2 {c\ a c ja +B..c.)+U J^n it n 4 +X J^(ct r c^ .+H.c.)(n<- ff +%_ <T ; 

(1) 

In ID, the model can display a phase with dominant singlet superconducting 



(SS) correlations, even for positive U [2l|, |26j, |27J, |28|. For X = t, the 



model has been solved exactly and there is a metal-insulator transition for 



increasing U 29, |30|, |3jJ. More recently, the role of entanglement in this 



quantum phase transition has been studied |32|, [33J. The response of the 
system to an applied magnetic flux indicates that the metallic phase is not 
superconducting [3l| . However, the ground state is highly degenerate in this 



phase at it is difficult to predict from the exact solution what happens when 
the degeneracy is lifted by a small but finite X — t. In any case, standard 



bosonization 20, 21 



26| and numerical studies 27] have provided a general 
physical picture of the behavior of the model, except at half filling. In this 
case, taking as usual only the leading terms in the lattice constant a, X 
disappears in the bosonization treatment [it enters as X cos(V7r), where v is 



the filling fraction [261. l21|]. Therefore, standard field theory predicts a SS 
phase for U < 0, and a spin-density wave (SDW) phase for U > 0, as in the 
usual Hubbard model. However, a charge insulator-metal transition driven 
by X at finite U c > has been found numerically [34| and later a quantum 
phase diagram has been derived which includes a BOW phase 35j . Recently, 



new numerical studies of the model at arbitrary filling identify regions of 
phase separation for X > 0.5t 



36 



1.2. Phase diagram at half filling 

The quantum phase diagram for v 
bination of different numerical techniques [35|, |36j. In Fig. [T] we repro 



1/2 has been obtained by a com- 



duce the phase diagram obtained by the method of topological transitions 



22. 37 



These transitions correspond to jumps in the charge and spin Berry 
phases which signal the corresponding transitions between the thermody- 
namic phases, and coincide with a corresponding crossing of excited levels, 
justified on the basis of conformal-field theory 37J, 13311 . T hese quantities are 
also related to charge and spin localization indicators 39_|, |4(J, |4l|, |42j used for 



example to characterize valence-bond-solid states in quantum spin chains 



3 




Figure 1: Phase diagram of the model at half filling obtained from the method of level 
crossings. Squares: charge transition. Solid circles: spin transition. The unit of energy 
is taken as t = 1. 
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431 ]. While changes in the Berry phase are proportional to changes in polar- 
ization, the spin Berry phase tensor provides a geometric characterization 
of the ferrotoroidic moment |42| . 

The critical values of U for the charge (U c ) and spin (U s ) transition have 
been calculated in systems of up to L = 14 lattice sites, and extrapolated 
to the thermodynamic limit using a parabola in 1/L 2 . Fig. [I] displays the 
extrapolated values for ^ X ^ 1. It is important to note that this interval 
can be extended to the whole real axis using symmetry properties of the 
Hamiltonian [27]. A change of phase of half of the sites [c\ a —> (— l) l cJ CT ], 
interchanges the signs of t and X, so that H(—t,—X,U) = H(t,X,U). 
Combining this with an electron-hole transformation (cJ CT — > Ci a ), one ob- 
tains at half filling 

H(t - 2X, -X, U) = H(t, X, U) = H(2X - t, X, U). (2) 

Then, if the critical ratio u(x) (U c /t or U s /t) with x = X/t, t > 0, is known 
in the interval ^ x < 1/2, it can be extended to negative values of X 
using the first Eq. ([2]). Similarly, the second Eq. (J2]) maps the interval 
1/2 < x ^ 1 onto x ^ 1. Explicitly 



if a: < 0, u{x) = (1 — 2x)u 



-x 



1 - 2x 



ifx ^> 1, u(x) = {2x — l)u 



X 



2x 



(3) 

It is interesting to note that the end points x — > ±oo are mapped onto 
x = 1/2. 

For U > U c (U < U s ), the system has a charge (spin) gap. For U < U c , 
the system is in the SS phase, while for U > U s , the system is in the 
SDW phase, according to the dominant correlation functions at large dis- 
tances. In between, for U c < U < U s , one has the fully gapped BOW 
phase. For small values of the interactions, the dominant correlation func- 
tions in each phase can be understood from field theory 



20, 21 



X > t/2, the SS phase displays incommensurate correlations [35 



For 
which 

can be qualitatively understood using a mean field approximation in one 



of the terms obtained from bosonization [35|], leading to a commensurate- 
incommensurate transition, with some similarities to the physics of the 
Hubbard model when a large next-nearest-neighbor hopping is added 44 



and some spin systems [45(. The spontaneously dimerized BOW phase has 
also been also found in the Hubbard model with alternative on-site ener- 
gies [37|, 46, 47, 48, 49, 50, 51, 52, 53], where it displays ferroelectricity 



371 . l5ll l52j . The presence of the BOW phase in this model was first pre- 



dicted using field theory and bosonization [46|, |47 



For X < t/2, the numerical results for the charge transition are consis- 
tent with U c = 0, as in the ordinary Hubbard model 35J. The accuracy of 
the results are not enough to establish if there is a kink or not at X = t/2, 
U c = 0. Arianna Montorsi has found that a good fit of the numerical results 
for t/2 < X < t is " 



54 



ifx ^ 1/2, u c = W2x - 1, (4) 

which is consistent with a kink at x — 1/2 = u c = 0, and has the nice 
property that when it is extended analytically to x ^ 1, it satisfies the 
second symmetry relation (131). The value U c = At for X = 1, is consistent 
with the exact solution 29|, |30|, |3l| To our knowledge, no justification of Eq. 
(jH) exists so far. 

It has been verified that the spin transition is of Kosterlitz-Thouless 
type [35j. In contrast to U c , the critical value for the spin transition U S (X) 
represented in Fig. [I] is smooth. For small values of X, it increases as X 2 , 
while for X ~ t/2 there is an inflection point. For X < t, U s > U c . While 
for X — t, U s = U c = 4t, the second Eq. ([3]) implies that for X > t, 
also U s > U c . Therefore, there is no crossing between U C (X) and U S (X) at 
X = t. The fact that there is a finite value of U s (l/2) and the second Eq. 
03j imply that U S (X) grows linearly with X for X — > +oo, in contrast to 
the \/X behavior for the charge transition predicted by Eq. 

1.3. Previous field-theoretical results 
As stated in Section 



II. H standard continuum limit field theory and 
bosonization fails at half filling because X disappears from the g>j coupling 
constants 



26, 211. In Ref. 35 



we have calculated vertex corrections to 
these gi using second order perturbation theory in the bond-charge inter- 
action X. The approach is similar to that done by Tsuchiizu and Furusaki 
for the Hubbard model extended with nearest-neighbor repulsion j56[, but 
for our Hamiltonian, Eq. ([1]) it is not necessary to introduce a low-energy 
cutoff. This approach led to the following critical U at the spin transition 



8A 2 



(5) 



n(t — X) 

This function lies below the numerical points in Fig. [H but seems to rep- 
resent correctly the limit X — > 0. However, unfortunately the prediction of 

6 



this approach for the charge transition is U c ~ U s /2 for small X, instead 
of U c = found numerically. In addition, with vertex corrections only, it 
is not possible to explain the nature of the incommensurate SS phase for 
X > t/2. 



In Ref. |35[ we have also considered in the bosonized theory, a term 
in next to leading order in the lattice parameter a, which couples charge 
and spin in a mean-field approximation. However, this approximation is 
questionable, and the quantitative agreement between the analytical and 
numerical results for the charge transition is poor. 

In this work we include all terms of next to leading order in a, and 
include them in a renormalization group (RG) treatment. This approach 



is superior to perturbation theory in X (as included in Ref. 35| through 
vertex corrections). Retaining all these terms leads to a lengthy algebra, 
but unfortunately selecting only a few of them, breaks the SU(2) symmetry 
and leads to wrong results. Since our approach is a weak coupling one, 
we restrict our study to X < t/2, which seems to be the more realistic 
regime of parameters. Our effort is rewarded by an excellent agreement 
with the numerical results for both critical values of U at the corresponding 
transitions. 



2. The field-theoretical approach 

2.1. The continuum limit 

In order to construct the low-energy field theory for the Hamiltonian 
Eq. (pQ), we suposse that both U and X are small. Therefore, in the Fourier 
development of fermion operators we retain only the modes near — kp and 
kp, where kp is the Fermi wave vector. Introducing a cutoff A << 1/a, 
where a is the lattice parameter, and calling L the lenght of the system, the 
local annihilation operator c na can be written as: 



E 



-ikpna 



Va [t 



E 

-A<fc+fcF<A 



i(k+kp)na < ikpna 



E 

-A<k-k F <A 



i(k-k F )na 
c u kcr 



no) + e ikpna 



(X 



na 



in the last step we have introduced the left and right fermionic fields ^) a -{x) 
and (x) respectively, by replacing the discrete lattice index n by a con- 

7 



tinuous variable x ~ na. This is possible because of the very small change 
undergone by sums in the second line of the previous equation, when one 
goes from site n to n + 1. 

Now we can undertake a gradient expansion for H by making the re- 
placement 

4>a±[x = (n+ l)a)] ->> 4j a ±(x = na) + ad x tp a ±(x)) (7) 
in all the terms of Eq. ([T]). For the hopping operator we obtain 

{c\ ia c n+la + c ] n+la c na ) ~ 2ia(-l) n (V^-Vv+ - ^ + VV-) + 

a 2 i[(ipl + d x ip a+ - d x ipl + ip a+ - ipl_d x ip a - + d x i\)\_^ a -) + 
(-1)" dM_^ a+ - iPUA-)} + 0(a 3 ), (8) 

where we have used kp = n /(2a). The number operator becomes 



n,. 



~ a p a+ {x) + p ff _(x)) 

+ (-1)" Ul + {x)^{x) + ^l_{x)^ + {x) 



(9) 



with p a+ = 4>i + ip a+ and p CT _ = Vi-Vv— 

By replacing by J — and taking into account that the integration 
of terms with an oscillating (— l) n prefactor vanish, one obtains for the 
Hubbard Hamiltonian [corresponding to the first two terms of Eq. (JTJ] the 
following form: 



H 



u 



iv F I dx I i/i + d x i/> ir+ - d x if>l + i/> a+ - if>l-d x i/> a - + d x ipi_4> a - \ + 



The first line corresponds to the usual free Dirac Hamiltonian with vf = at 
the bare Fermi velocity (later we shall use Vp = a(t — X), the Hartree-Fock 
value [21]) The terms with prefactors gi±, #2±, <?3J_, Qa± correspond to back- 
ward, forward two branch, Umklapp and forward one branch, respectively. 
While all these constants are equal to aU, they might run independently 
under a renormalization group (RG) flow. Moreover, in units in which 

8 



(10) 



h = 1 and Vp = 1, the couplings are dimensionless. Now, if a coupling 
constant gi has units E d ~ Ai , Aj is known as the scaling dimension of the 
corres pon ding operator Oj, where d is the spacetime dimension (2 in our 
case) [57|]. Therefore all interactions in Eq. f JT0|) have scaling dimension 
2, they are marginal operators. It is known that depending on the sign of 
the Qi which correspond to the charge or spin sector of the theory, they can 
become marginally relevant or irrelevant. The first case leads to a charge or 
spin gap [58( (see also Section |2"^|) . Let us see how this situation is modified 
by inclusion of the correlated hopping term [the last one in Eq. (1T1)] 

The sum of the number operators in (pQ) has the following gradient ex- 
pansion: 

(n nW + n i+m ) ~ 2a(p^+ + Pa-) + a 2 [d x (pa+ + Pa-) - 

(-l) n (dM-^ + + 4 + ^-))}+0(a 3 ) (11) 

multiplying (jHj) by f TTTj) we see that the terms quadratic in a are oscillat- 
ing and vanish under integration. This is the result anticipated in Section 
11.11 This means that no scattering at the Fermi level is generated by the 
correlated hopping interaction. We should include term up to 0(a 2 ) in the 
Hamiltonian. We obtain: 

H x = i dx^2 {#4± [-^4+^+4+^+ + 9 b V'J-'0<t-V4-VV- - H.c] 

^ a 

g' 2± [-dxi>t + ^a+^l-^a- + d x ipl_ip a -ipl + ip w+ - H.c] 
g' 11 \- d x^l+^a-^l-i^+ + <9^_Vv+V4+^- - H - c -] 
g' 3 A-d x il>t-^+4-th+ + 0,^+^-4+^- - H.c.]} (12) 

where all g\ = a 2 X. The essential differences between the field theory given 
by Eq. ( I12p and the one given by Eq. fflUj) is the non local nature of 
the interaction arising from the derivatives. In /c-space this corresponds to 
scattering of electrons which are near but not on the Fermi surface. 

Note that the coupling constants g\ has dimension of the inverse of en- 
ergy. Therefore each term of Eq. (fl2|) has dimension 3, they are irrelevant. 
As only these irrelevant operators appear in the low energy fermionic field 
theory of the correlated hopping term Hx, one might be tempted to con- 
clude that there is no contribution of these terms to the physical behavior 

of the system, in contrast to the numerical results discussed in Section 11.21 

9 



How could we account of this situation with our field theoretical analysis? 
In fact, from a renormalization group (RG) point of view, all the operators 
allowed by symmetry should be included in the effective theory. The fact 
that operators present in Eq. (11 01) were not obtained in the derivation of Eq. 
(j!2p means that in the initial conditions, the different g do not depend on 
X. However, they could acquire an X dependence by the couplings between 
g and g' when the RG flow evolves. 

We can take advatage of the well known studies of thermal critical phe- 



nomena with RG [59| to further understand this issue. For this case we 
know that the irrelevance of an operator means that the critical exponents 
are not affected by its presence in the Hamiltonian. However the critical 
temperature does depend on this operator. In our case the presence of ir- 
relevant operators will be crucial to determine the boundaries in parameter 
space of the different phases, where the spin or charge gap opens. 

Finally, we note that the SU(2) invariance of the ordinary Hubbard 



model implies that under RG flow, g2± = g±± remains [58|. Similarly, it 



is shown in Appendix A that g' 2± = g' 1± is required to keep the SU(2) 



invariance of the full Hamiltonian. 
2.2. Bosonization 

Bosonization is a powerful technique to analyze interacting one-dimensional 



fermionic systems |58|]. Some of the interacting terms in the fermionic 
Hamiltonian become free non-interacting terms in the bosonic Hamiltonian. 
The remaining terms contain in general cosines of the bosonic fields. Their 
effect can be studied by a perturbative implementation of the RG method. 
If in the RG flow the coefficient of a cosine decreases, the fixed point corre- 
sponds to a trivial theory of free bosons with known properties. When the 
RG flow goes to strong coupling the coefficient of a cosine increases. The 
fields are trapped in a minimum of the free energy and the different phases 
can be characterized by calculating the classical value at this minimum of 
the bosonic operators corresponding to the physical observables. In our 
case, the RG analysis is more involved, but as we shall show, it leads to a 
tractable theory and correct results. 

Let us therefore resort to a bosonic representation of the fermionic theory 
of Eq. (j!2p We use the following bosonization formula for the left (— ) and 
right (+) fermions 6(| : 



^ a± (x) = ^=:e^: (13) 
10 



V2 



- e =F<*„± ( 14 ) 



Equation ( 1T3|) is normal ordered and therefore does not contain a somewhat 
uncomfortable short range cutoff a, F is the Klein factor and L the length 
of the chain. Eq. (114"]) arises from Eq. (fl3|) by the explicit expansion of CT ± 
in term of the boson creations ib a ±) and annihilations (b\±). It is given by 



58|,l60 



_^ ^fipx—ap/2 ^ g±ipx— ap/2 

j±(x) = TiY] =— bl±±iY] &<x±, (15) 



<4± ^ ct± 

where (<£> CT ±) are the creation (annihilation) part of the field <p a ± and 
P = l^T - ^ e Produce the charge and spin bosonic fields pc and 4> ps 

(p = +r~) 

, _ (0gt + 0p+) , _ (0Pt ~ 0pJ.) 

<P P c~ ^ ,<P P s- . 

We also introduce phase fields (f> m and 9 m {m =f, J,, c, s) 

ji 0m+ + 0m- fl fym- ~ 4 > m+ n r?\ 

v m = [il) 



(16) 



2 2 



The line before the last in Eq. (fT2|) bosonizes as: 

9l±) (- d ^l+^a-tpl-lh+ + d x lfjt_^ a+ 1pl + ^ w - - H.C.) = 

^[e^a^+e" 2 ^ + e-^d^-e 21 ^} - H.c. 



2 



(2vrct) 

4\/2z 4\/2z 

cos(2V20 s )c' a; c = W —: cos(2v / 20 s ): (18) 



(27ra) 2 v v ™ xrc L 2 

In the last line we have normal ordered the cosine using Eqs. (j!3p and ( 114)) . 
We also have: 

^ cos(2v / 20 c )«9 :c c = ^:cos(2V20 c ):^0 c (19) 
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The bosonization of g' 2 j_ and g' A± terms is a little more subtle. It is convenient 
to come back to the lattice version of the derivate with respect to 
given in Eq. (17j). We have: 



%l) rT+ %_ijja- + 



+ ( V>£- (a; + a)^- (x) - (x + a) J 

We use Eq. ([13]) to bosonize the first term into: 



(20) 



j_ 

La 



Va+(x)^pa+(x + a) - ^a+(x + Cl)ljj a +{x) ) = — 



e i(j> a+ (x). . e -i0 CT+ (x+a). _JJ c 



2,„t 



Li \ 



H.c. 



2na 2 



a J \2 



i(l - iad x (pl + - zy<9^ + - ^-{d x ipi + 
x [ \ - iad x ip a+ - i—dlLp u+ - — (d x Lp a+ ) 2 ) - H.c 
2(l - a 2 d x (pl + d x ip a ^ 



2na 2 



l + y:(40 CT+ ) 2 : 



(21) 



In the second equality we have used the identity 60j e A e B = e B e A e^ A,B ^ 
being the commutator a c number which could be calculated by the explicit 
expansion given in Eq. ( fT5|) with a = 0. The commutator becomes: 

gipx g—ip' (x+a) 



[ Va+ {x), V l + {x + a)} = 7777 

p,p'>o V p p 



-1 pa 



S P, P > 

2na. 



Then, 



J^<T+(xO,^L(x+a)] 



iL 



2rra ' 



(23) 



the value previously used. 

Proceeding in a similar way one finds (to be used later) 



,[tp a -(x),tpl_(x±a)] _ _|_ 



iL 



2na 



(24) 



The term into the parenthesis in the last line of fl20l) can be bosonized 
by similar steps. The result is: 



- ifjl_(x+a)ilj a _(x)-iljl_(x)ilj (T ^(x+a) 



-1+-:(^0 CT _) 2 : 



a \ J 7rcr 

Taking into account that the normal ordered densities bozonize as: 



(25) 



1 



: p a ±: =: V£±Vv±: = ~7r : d x4>a±- 

Z7T 



we obtain the bosonized expression of Eq. f l20 



2(vra) 



2(tto) 2 

V2~(d x (fr c+ + d x (p c ^j 



(d x (f) a+ + d x (j> a -) - y ( (d^+fd^- + (d x (f) a -) 2 d x (f)^ 



^ ({d x (j) c _) 2 d x (j) c+ + d x (j) c+ ((d x (j) s J) 2 - 2d x (p s „d x (p s+ ^ + 
d x (j) c - ((d x (f) c+ ) 2 - 2d x (p s -d x (f) s+ + (d x (p s+ 

Quite similar steps leads to the bosonization of the g' A term. We find: 



(26) 



2(7ra) 2 



(d x </> a+ + d x (f> a . 



(d x 4> a+ ) 2 d x fe + + (d x <p^) 2 d x (h- 



2(nay 



V2~(d x (p c+ + d x (f) c ^j 



2V2 



{d x <p c -f + {d x <p c+ f - d x <PUdAs~? - <ho, .{i) AK .Y 
13 



(27) 



Collecting the different pieces, going to an imaginary time r = it and defin- 
ing complex space-time coordinates (z = vpT + ix,z = VpT — ix), where 
v f — o^t — X) is the Fermi velocity (starting from a Hartree-Fock decou- 
pling [21]) we obtain the following expression for the part of the action 
proportional to X. 



S x = a(G\J d 2 rO[{r)+G' 3 j d 2 rO' 3 d 2 r + G' 2 J d 2 rO' 2 d 2 r + G' i J d 2 rO'^d 2 r) 

(28) 

where G' a = g' ±a /(aiTVF) and d 2 r = vpdxdr. The different operators in Eq. 



28J) are: 



0[ 



O' 



OL 



O' 



i — ■ cos(2v / 20 s ): (<9 Z C _ - c^0 c+ ) 



— ^- 



L 2 
.2tt\/2 

2na 2 



: cos(2\/20 c )(^0 c _ - <3^> c+ ): 



(dz4> c+ - d z (f> c S) H -=- ( : (d z (f) c -) 2 dz<f) c+ : 

v / 4V27T V s „ ' 



Oil 

+ ck^cV- (d z (f) s ^) 2 : + 2& t -(j) c+ d z (f) s -d z -(f) s+ 



O' 



O' 



o' 2A 



( Qz<\>c-- (dz<Pc+) 2 -- + 2d z (p c ^d z <p s ^d Y (p s+ + d z <p c ^. (<9^> s+ ) 2 : ) J 
V s v ' v V ' v v J 

Or, K Oo ft Oiy ? 



^2_ 

2a% 

(<9 2 c _) 3 - (^0 C+ ) 3 - <9 2 <M<9 2 S _) 2 + ^0 c+ (^0 s+ ) : 



(^0 C + - d z 



AV2-K 

To obtain Eq. ( 128]) we have: 



(29) 



1 . Included the normal order of each bosonic operator assuming that the 
original fermionic operators were already normal ordered. This is a 
prerequisite for the bosonization to work (soj ] . 

2. Taken into account that d x = i(d z — dz) and 

3. that the right and left bosons depend on z and z respectively. I.e. 

4> m +{z) and (f) m -{z), (m = c or s). 

14 



This last fact arises when the explicit time dependence of the bosonic 
creation and annihilation operator is deduced from the Heisenberg equations 
of motion using the free bosonic Hamiltonian H = vp J2k>o ^Knk+^mk+ + 
v f J2k>o kb ] mk _b mk „. One obtains & mfc ±(r) = exp~ VFTk b mk ±(0), and b ] mk± (T) - 
eX pV F rk 6^ fc± (0). Plugging these expressions in equations like ffT5l) one ob- 
tains: 



e —ipx + pvft e %px — pvft 



t m + 

— pz pz 



e ~tpX — pVpT e ipX + pVpT 

- % H ^= b mp _+tJ2 -= b j mp _ (30) 

n p >0 V P n p >0 V"V 



S 



<Pm- J 



The total action is S = Sh + Sx, where Sh is the usual bosonized version 
of the Hubbard model of Eq. ffTOj) 



Sh = 2~ / d z<f>c-dz(t>c+d 2 r + — / d z (f) s ^d^f) s+ d 2 r + 



G?i y rf 2 rOi(r) + G 3 / d 2 r0 3 d 2 r + G 2c J d 2 r0 2c d 2 r + G 2s J d 2 r0 2s d 2 r + 
5v s [ d 2 rO is d 2 r + 5v c [ d 2 rO ic d 2 r (31) 



with 



01 = 


-^:cos(2\/20 s ): 


(32) 


o 3 = 


^:cos(2V20 c ): 


(33) 


o 2c = 


47T 


(34) 


o 2s = 


-— d z (f) s -ck4>s+ 

47T 


(35) 


o 4c = 


^((^ c+ ) 2 +(^ c _) 2 ) 


(36) 


o 4s = 


-^((^0 s+ ) 2 + (^0 s _) 2 ) 


(37) 



Sv c (Sv s ) renormalize the charge (spin) velocity. Operators 02 C arid O^s 
(and 4c and O is ) appear together in the bosonized theory of the Hubbard 
model Eq. (imp , where only interaction between electron of different spin 
are taken into account. They are independent operators in the general case 
where interaction between electron of the same spin are included. For the 
Hubbard model, the values of all couplings are G a = Uafiirvp)- 

2.3. The renormalization group equations 



Following Ref. [6jJ, the RG equations for the coupling constants r a (G 



or G' a ) present in the action is 

^ = (d - A 7 )r 7 - Yl °yv>. (38) 

a/3 

where d — 2 is the spacetime dimension of the system, A 7 is the scaling 
dimension of the operator related with T 7 , Sd is the area of a sphere of unit 
radius in d dimensions (27r in our case), the CZg are the coefficients of the 
following short-distance Operator Product Expansion (OPE): 

q / r+r' \ 

0' a (r)0'p{r') = J2 G *P ) r _ r > \a!+a p -A-, + more irrelevant operators, 

7 

(39) 

and A is a number of order one, which comes from our definition of the 
short distance cutoff as a/ A (the exponent of A in Eq. (|3"8"]1 comes from the 
integral of Eq. fl39|) with respect to | r — r' \ in d space-time dimensions). 
The OPE's between two O a operators are already known from the RG 



equations of the ordinary Hubbard model [58[]. The OPE's between one O a 
and one O'p operator give another 7 and have a prefactor G a G'p ~ UX. 
We note that for small X, this product is of order X 3 on the spin transition 
and of higher order or negligible on the charge transition. We have neglected 
these OPE's. This is partially justified by that fact that they generate 7 
operators which are irrelevant, while as we show below the OPE's between 
two 0' a operators generate marginal O a operators. A deeper justification 
in given on symmetry grounds: expressing the first Eq. ([2]) in terms of the 
Hartree-Fock hopping t = t—X [which is invariant under the transformation 
ctj — y (— iyc ia ] one has 

H(l -X, U) = H(t, X, U). (40) 
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This means that for small X, there can be no terms of order UX in the action 
which correct the Hartree-Fock results. Therefore, the generated operators 
in the OPE's between one O and one 0' a should introduce corrections of 
higher order. 

Now let us discuss the different operators that could arise from the 
OPE's between two 0' a operators. There are some cases where these OPE's 
give operators of dimension 4 or higher. This is for example the case of the 
OPE between 2 3 and 27 or in general between two operators included 
in O2 which contain less than two fields in common. There are other cases 
where the denominator does not depend only on the distance between the 
two points under consideration but have factors of the form (z' — z)~ 2 + 
(z' — z)~ 2 . This gives rise to a periodic function in the relative angle and the 
integral in the angular part of r — r' [which was performed to arrive at Eq. 
(1551) ] vanishes. This is the cases of the OPE between 0[ and (0' 2A + 0' 26 ). 
Finally there are cases which produce operators of the form (j34p and (155]) . 
They simply renormalize the charge or spin velocity and will not be taken 
into account in our treatment The remaining OPE's are displayed in the 



Appendix B| From this appendix we have: 



U 2'V 


- r 1 


1 


r 3 

°2'3' 


- r 3 


1 


si2s 
Ly 2 , 2' 


r~i2c 
— °2'2' 


1 

71 









(41) 

From Eqs. ( 1551) and ( 14T|) . the usual RG equations for the Hubbard model 
become modified as follows. For the charge sector 

^ = G\- \\G' 2 )\ ^ = G 2c G 3 + A 2 G' 2 G' 3 , (42) 
for the spin sector 

^ — ~G\ — A 2 (G" 2 ) 2 , -— - = — G 2s Gi — X 2 G / 2 G / 1 , (43) 
and in addition 

fir" 

— ^ = —G' . (44) 
dl a K ' 
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2.4- Analysis of the RG equations 

Taking into account the initial conditions, Eq. (144 j) can be integrated 
immediately giving 

C = — «p(-Q. (45) 

TTVf 

Replacing this equation in Eqs. ( I4"2j) . one obtains two coupled differential 
equations for the charge sector 



^ = Gl - Ae~-\ ^ = G 2c G 3 + Ae^\ with A = (^) ' , (46) 
dl dl \ ttvf J 

with the initial conditions G2 C i} = 0) = G 3 (l = 0) = Ua/ (ttvf). 

It is known that for A = 0, the flux continues along the separatrix 
Gic = G3, and goes to infinite G{ (charge gap) if U > 0, and to G 2c = G3 = 
(gapless case) if U < 0. While an analytical solution for A / seems not 
possible, it is clear that the effect of A is to push the flux perpendicularly 
to the separatrix, favoring larger G 3 and smaller G 2c - This does not modify 
the final result that the critical value of U which separates the regions of 
diverging or vanishing Gs(l — > +00) is U c = 0. We have confirmed this by 
a numerical study of Eqs. (H6|) . However, as a difference with the Hubbard 
model for which the flux is on the separatrix, in our case, for U < U c 
(when G 3 flows to zero), G 2c converges to a negative value. This leads to a 
correlation exponent [58| K c ~ 1 — G 2c larger than 1. As a consequence, the 
singlet superconducting (SS) correlation functions, which decay as d~ l l Kc at 
large distance d dominate over the charge density wave (CDW) ones, which 
decay as d~ Kc 2o|, [HI, [HI] . In the Hubbard model, for U < 0, K c — 1 and 



both SS and CDW correlations decay as 1/d. 
For the spin sector, the RG equations become 

^ = _ G? _, e - a ,^ = _ GlG2 ,_, e - a „ (47) 

with the initial conditions G2 S {0) = G*i(0) = Ua/^Vp). 

It is clear that the flux of the RG equations remains on the separatrix 
G\(l) = G2s(l)- Therefore, both equations fj4T|) reduce to the same equation 
for Gi = G*2s = G. Changing variable z = \/Ae~ l , this equation takes the 
form 

dz lg 



Its solution is given in terms of Bessel functions 



where the constant C is determined by the initial condition G(y/A) = G mi , 
giving 

-G ini J (VA) + VAJi(VA)' 
Mathematically, for I — > oo (z — y 0), Eq. fj49|) converges to zero. How- 
ever, it may happen that G(zd) diverges for some intermediate value Zd 
(0 < z < \fA)i jumping from — oo to +00 as z decreases. This means phys- 
ically that at an intermediate scale determined by z^, the solution flowed 
to the strong coupling fixed point at which a spin gap opens. The limiting 
value of Zd for which such a behavior takes place corresponds to Z4 — > 0. 
Since for small values of the argument Y (z) ~ (2/7r) ln(z) and Jq{z) ~ 1, a 
diverging G(zd) for Zd — > implies a zero in the denominator of Eq. (T49"|) . 
and the initial conditions should be such that C also diverges in this special 
case. 

For small values of A (as we have assumed in our whole treatment), 
there is no divergence in G(z) if C is negative. From this reasoning and Eq. 
f l5T)|) . we obtain the following condition for the opening of a spin gap: 

M</A) 2 16' ^ ' 

where the last member was obtained from a series expansion of the Bessel 
functions. 

From Eq. (IBT]) and using G ini = Ua/invp), we obtain the following 
critical value of U for the opening of the spin gap 

VAMs/A) 

u - = ™ F ^7ir (52) 

or approximately 

The final path taken by the RG flow in each sector, determine the na- 
ture of the resulting phases. As discussed above, for U < U c , SS correlations 
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dominate. For U > U s spin-spin correlations are the largest at large dis- 



20, 21, 58 



tances as in the usual Hubbard model 

The phase in between, for U c < U < U s is characterized by the presence 
of both gaps, and the RG flow in each sector leads to G 3 — > +00 and G\ — > 
—00. To minimize the respective cosine terms in the action [See Eqs. ( IHTj) . 
fl32l) and (133]) ]. the fields are frozen at the values 2y/2(p c = n and 2y/2(p s = 0. 
As a consequence, the system has a spontaneously dimerized bond-ordering- 
wave (BOW) phase with long r ang e order. The order parameter which takes 



a finite value on this phase is [20j, |21 



Obow = ^(- 1 ) i (4+i jCT ^ CT + H.c.) ~ sin(v / 20 c ) cos(v / 20 s ). (54) 

3. Comparison with the numerical results 

The field theoretical result for the charge transition U c — obtained 
in the previous Section, agrees with the numerical results, presented in 
Section 11.21 As explained in Section ll.3[ this result is not obtained if the 
initial values of the couplings of the Hubbard model ((j2c(0) and Gs(0)) are 
corrected by vertex corrections in second order in X before bosonizing. We 
do not have a physical explanation for this. 

To compare the critical value of U for the spin transition, we take A = 4, 
which using Eq. fl53l) and recalling that vf = a(t — X), leads to the same 
result as that obtained from vertex corrections Eq. (j5J) for small X. This 
leads to A = {4X/[n(t — X)]} 2 . Replacing this result in Eq. (1521 . we obtain 
the function U S (X) that is represented in Fig. [2J The agreement with the 
numerical result up to X/t ~ 0.4 is excellent. 

4. Summary and discussion 

We have studied a field theory for the Hubbard model with small bond- 
charge interaction X at half filling. While usually, it is enough to consider 
in the action only terms linear in the lattice parameter a, in our case it is 
necessary to include terms of order a 2 to obtain meaningful results at half 
filling. These terms can be classified in a similar way as the linear ones in 
terms of different processes in a "g-ology" treatment (forward one branch, 
forward two branch, backward and Umklapp) but contain derivatives of the 
fields in the space direction. 
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Figure 2: Solid circles (squares): critical value of the spin (charge) transition U s (U c ) as 
a function of X obtained from the method of level crossings. Full line: solution of the 
RG equations given by Eq. (|52l) . The dash dotted line at U c — signal the boundary 
between the gapped and gapless charge phase as given by Eq. (|46|) . The unit of energy 
is taken as t = 1. 
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We have obtained the RG equations of the different couplings using 
Operator Product Expansions. While the treatment of the new terms is 
awkward, most of them should be retained to keep the spin SU(2) invariance 
of the model. 

According to the dominant correlations at large distances, the phases 
of the model can be classified as a singlet superconducting (SS) one for 
U < U c , a bond ordering wave (BOW) for U c < U <U S and a spin density 
wave (SDW) for U > U s . The boundaries between the phases correspond 
to a charge transition for on site repulsion U = U c and a spin transition at 
U = U s , which correspond respectively to the opening of a charge gap and a 
closing of the spin gap as U increases. For the former transition we obtain 



U c = in agreement with previous numerical studies for X < 0.5 [35] . With 



only one adjustable parameter, we also obtain a very good agreement with 
the numerical results for U s if X < OAt. 

To explain accurately the dependence of U c for X > 0.5, it is necessary 
to go beyond our approach, possibly including more irrelevant operators. 

As stated in Section 11.11 the model is an effective one-band model for 



a variety of physical systems, in particular optical lattices [121 . 1 131 . [14], [15 
In these systems, U can be varied over the whole range, including its sign, 
through tuning of the external magnetic field B. It is also possible to change 
X/t by 20%. Therefore, adjusting the filling to one particle per site, it seems 
in principle possible to tune the parameters in such a way the ground state 
of the system is in any of the three phases: SS, BOW or SDW. 
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Appendix A. SU(2) spin rotational symmetry 

In this appendix, we derive the relations among the constants g' which 
leave the theory invariant under spin SU(2) transformations. 

Under an SU(2) rotation each spinor i/; pa , (p = +, — ) transform as ip'^, = 
U a i a ip a (the sum over repeated indices is implicitly understood). U is an 
SU(2) matrix given by: 
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u =U (A1) 

a and h are complex number satisfying |a| 2 + |6| 2 =l. We require the 
invariance of the g-ology Hamiltonian (IT2"1) under this rotation. We start by 
undertaking a rotation of the sum of the first terms of the second and the 
third line of (IT2|) . The transformed expression is: 

~ K,a U cra 2 U*a U aa 4 ( d x 1pL + 1p a2+ 1p*1p a4 - - ^4+^4+4-^2- + 'WaV'ai+V'o 



Vio- < - / o-CT 2 < - y ( T3CT w O-CT4 I ^1^(71 + ^2+ ro- 3 - rCT4- ^ *r (Ti + t 0"4+ 7(73 — r 0"2 — ~ u cr 2 (T3 r O"! + ro-4- 



+ ^CU^^A+Vw (A.2) 




In the second term of the last line we have interchanged the dummy indices 
02 with <7 4 . Now it is possible to evaluate the difference of the products of 
U matrices involved in flA.2j) . We have: 

fi(T2V4.(.U aiU2 U asW2 — U ai - d : 2 U IJ3a2 ) = 5 a2 a 4 5 a . i a 1 (8 Uia2 — 5 a2U3 ) (A. 3) 
Inserting in (I A . 2 [) we obtain: 

-(3,4+^+4-^- + 0*4+^-4-^+) (a.4) 

This is the sum of the first terms with the coefficients and g' 2J _. Therefore 
we have shown that the sum of these terms is invariant under an SU(2) 
rotation. With the same procedure we can show indeed that the sum of 
all the terms multiplying g' 1± and the ones multiplying g' 2± transform into 
themselves by an SU(2) rotation. This implies that under the condition 
g' 1± = g 2J _, the Hamiltonian Hx of Eq. ffl2l remains invariant. 

Appendix B. OPE's between O' operators 

In this appendix we give the expressions for all non-vanishing OPE's 
between any two O' operators. The first one is: 

0' 2 . 7 (z,z).0[(z',7) = -i=-:^(^)(^0 s+ (^)) 2 : 
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x ^?:cos(2v / 20 s (z , ,^)(^'0e-(^) - d*<j) c+ {?)): = 

^2 (: cos(2V20.(z', z'))(^'^+) 2 : + =rr= : sin(2v / 20 s+ (2 , ,^))^. a 



+ ... (B.l) 



1 \ : cos 



.L 2 J \z'-z\ i "' 2tt\z'-z, 

To obtain the previous result we have first normal ordered each product of 
operators containing fields at different points. Then we have expanded the 
resulting expressions for z' near z. Let us show how we have proceeded step 
by step : 

d z <f> c 4z) d z <f) c 4z') = (d z <pUz)+d z <p c 4z))(d z ,d-(z') + d z >iPc-(z')) = 

d z (p c ^(z)d z/ ip c ^(z') + d z tp\_(z)d z iLpi_(z') + d z <pl_(z)d z ,(p c -(z') + 

dj<pi_tf)d z <p^(z) + d z dA<P^(z),<pLV)] (B-2) 

and: 

[ ( Pc-(z),(pt_(z')} = 
„—pz„p'z' „^f-n v (z' —z) 

^ = ^> -log(l-e-<"'>) 

v s p,p' |e ^ 1 j |< 1 

v «, ' 

t'<T 

-log(^(>- z')) 

L»|*'-*| 

a,a,[^(*),rf_CO] = - {z ^ z , )2 ( B - 3 ) 

The third equality has been taken from Ref. [62]. The condition r' < r 
implies that the operator should be time ordered in decreasing order from 
left to right. Including Eq. (1B.3j) in Eq. f lB.2j) we obtain: 

<9 2 c _(z) d,<f>^)=:d s (i>^z) d z ^{z'):- . 1 (B.4) 

(2; — z J z 

which is the result used in Eq. (IB. ip . Regarding the normal ordering of the 
product : (d z (f) s+ (z)) 2 : : cos(2\/2(f) s (z' , z'):, we use the following basic OPE: 
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x ^d^ s+ (z) + iX [^ s+ (z), <pl@) } e iX ^ (B.5) 



1 



2(z'-z) 



where we have used Ae B = e B A + [A, B]e B and the commutator was calcu- 
lated as in Eq. (1B.3[) . Therefore, we have: 



dzip s +{z): cos(2V2(f) s (z' ,z')): 



: cos 



(2>/2>.(z',*'))3*¥>,+(z) : 



V2 



\z - z 



-:sin(2v / 20 s (^ , ,2 / )) : ( B -6) 



(the last sign changes if one permutes sin and cos) and 
(^ + (^)) 2 :cos(2v / 20 s (^ / ,^)): = 

y/2 



d^p s+ {z)(: cos(2\/2> s (z' ',*)): dg<p a+ (z) - ,_, : sin(2v^0 8 (^,^)): 
V ( z -z) 

(:cos(2\/20 s (^,^)):^+(^) - ( -/^ : sin(2v^0 a (^,^)): )d g <p a +(z) 
-7=^=r (: sin(2v / 20 s (z', z')): d^p s+ + 



V2 



(z' - z) 



cos 



(2V2<p s (z> \z')): 



:cos(2V20 s (/,z / )): ((«%^)) 2 - ( _/_ )2 ) 



2\/2 

^ — — : sin(2v / 20 s (2; / , Qtys+i 



(B.7) 



The required OPE is: 

:(dz(P s+ (z)) 2 : :cos(2V20 s (z',I / ): = 
(<3^ s+ (z)) 2 + 2(<3^ + (z))(«3^ s+ (z)) + (<3^! + ) 2 ] : cos^v^z',^): = 

:cos(2v / 20s(^ , ,^ / ))(^'0 s +) 2 : 

n (B.8) 



2\^2 2 

" : sin(2A/20 s+ (z',2 / )): -— — — : cos(2\/20 s (y, 



z' — z 



\z - z 



which is the other OPE used in Eq. fIB.ip . 
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An analogous reasoning leads to the remaining non-vanishing OPE's: 



O' 2 . z {z,z)0 x {z\z<) = -^—■.d, ( f )c+ (z)(d z( f )s 4z)) 2 : 

4V27T 

x— —: cos(2V2(p s (z' ,z')(d z ,(i) c -(z') - d^>(f) c+ (z)): = 

2^2 



2 ^ (: cos(2v^0 s (/,^))(^0s-) 2 : sin(2v / 20 s+ (/,^))^0,-: 



+ ^ z , z z y cos(2V2(f) s (z' , z'))\ ) ( + ^_ + ^0 C+ ^0 C _-: ^0 c+ ^/0 c+ : ) ~ 

/ 1 ycos(2x/20 s (^) 1 Oi , (nQ , 

0^{z,?)O z WJ) = -^-:^0 C+ (^)(^0 C _(^)) 2 : 
x^?: cos(2v / 20 c (^ / ,^))(^0 c -(^) - M^ffl): ~ 

-feJ |^7p — - = "27R^7F + - (R10) 

Oi 5 (z,z)(y 3 ^^) = --^ r :d z <t> c 4z)(d z - ( f> c+ (z)) 2 : 

4V27T 

x^?: cos(2v / 20 c (^ ! ^ / ))(^'0 c -(^) - Wc+CO): = 
-(^)^c-(:cos(2v^0 c (^O)(^ 

- 2 : cos(2v / 20 c (^, ^))dAc-- ) - 
(z' — zy ' 

-feJ — |V^7F — - = "2^R^7F + - (R11) 

OUz,z)O f 2A (z',z') = -^-^:d^ c+ (z)(d z <j )s 4z)) 2 : 



x ^0 c+ (^)9 z /0 s _(2; / )^0 s+ (2; / ): ) = 



^(:(^ W )»^ (J 0:-^=W) 
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1 dz<j> a+ {z)d z (j) a -{z)_ + 1 2s | ^ B 12 ^ 



8tt 2 I z - | 4 

1 / 1 



2 . 7 (^)0 2 . 6 (z',z') = --(^^(^(S^+C*)) 2 :) 
x ( — -: dAc-{A Q As-{A^s^)'- 

v V 27T 



167T 2 V r*^~„ r,TV~ /■ ( - _ -, )2 

x(:<9 z c „(2:) d z ,(j) c -{2!):- -j-^ )d z /0 s _(z') = 
1 3^ a+ (z)0 z a _(z) 1 2s 



8tt 2 l^-^'l 4 " 27r|z'-z i 

_ _ 1/2 _ 

2 . 4 (^^)C ) 2 . 6 (^,^) = — \^-^:d z -(f) c+ (z)d z (f) s -(z)dz(f) s+ ( 



x 



7T 



: ^0 c _(z')^/0 s _(z / )9j'0 s +(^) : ) - 



1 ^0 c+ (z)9 z c _(2;) 1 2c 

+ - = ^Z|-3 Tu+- ( B - 14 ) 



8tt 2 U-z'I 4 " 2tt\z'-z 



Oi 2 (z,z)0' 2 _ 5 (z',z') = {d z <P c ^z)f<k<P c+ {z):) (:^0 c _(/)(^0 c+ (^)) 2 :) = 



32tt 2 V v ^ c_v /y ~*™-v~'- ( Z ,_ Z ) 2A - c ^ v yv ^ Ci " v " • 

1 <9z0e+(X>d Z 0c-(^) 1 Q2e , R 1 v 

U_ 2 ' |4 + - - -.4 + - 
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